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A Subvortex Technique for the Cloee Approach 
to a Dlscrctlaed Vortex Sheet 

Brian Maskew* 

KASA Asiee Research Center, Moffett Field, Calif. 9A035 

The closed-approach problem associated with flow calculation methods 
based on vortex^lattlce theory Is examined nm^rlcally using two-dliMnslonal 
discretized vortex sheets. The analysis first yields a iiear-fleld radius of 
approximately the distance apart of the vortices in the lattice; only within 
this distance from the sheet are the errors arising from the discretization 
significant. Various modifications to the discrete vortices are then con- 
sidered with the objective of reducing the errors. This study leads to a 
near-fleld model In which a vortex splits into an increasing number of sub- 
vortlces as It Is approached. The subvortices, whose strengths vary linearly 
from the vortex position, are evenly distributed along an Interpolated curve 
passing through the basic vortices. This subvortex technique can be extended 
to the three-dimensional case and Is efficient because the number of vortices 
Is effectively Increased, but only where and when needed. 


The author wishes to thank Mrs. Opal J. Lemmer for her assistance with 
the computer program for this work. He Is also grateful to Dr. V. J. Rossow 
and Messrs. R. T. Medan and V. R. Corslglla of the Large-Scale Aerodynamics 
Branch at NASA Ames Research Center; the paper has undoubtedly benefited from 
their Interested discussions during the course of this work and also from 
their comments on the first draft. 

*NRC Research i^soclate on leave of absence from Hawker Siddeley Aviation 
Ltd., (Brough), North Hxmberslde, England. 


- 1 - 



1. Introduction 


Velocity distributions calculated close to a discretized vortex sheet 
show distortions because of the singular behavior of the Induced velocity 
field near each vortex. The problem is associated with all vortex-latt ice- 
based methods (e.g. » Refs. 1-5), but it has been circumvented in the past by 
calculating near-field velocities only at special points, e.g., midway between 
the vortices, or by Interpolation between velocities at nearby "good” points. 
For calculations Involving multiple vortex sheets,'**® the near-field problem 
often requires that adjacent lattices be made to correspond across the gap 
between the sheets. However, such a solution is not practical for recent 
developments in vortex-lattice methods which incorporate iterative calculation 
schemes for trailing vortex relaxation®*® (i.e., force-free wake). Although 
the new methods have proved to be very versatile in general, close-approach 
situations involving multiple discretized votex sheets require careful treat- 
ment, and Ideally, the near- field problem should be removed. 

The purpose of this Investigation, therefore, is to develop a near-field 
modification for the discrete vortices which will allow velocities to be cal- 
culated anywhere in the flow field irrespective of the proximity of discre- 
tized vortex sheets. Such a modification would enhance the versatility of 
vortex- lattice-based methods. 

In the Investigation, a near-field region is first defined and its 
extent examined numerically (Sec. 2) before considering two near-field models 
(Sec. 3) fdilch lead to the subvortex technique. The latter is described in 
Sec. 4, and an application is given in Sec. 5. The analysis is perforrod 
using two-dimensional situations, but the application to the three-dimensional 
case is taken into consideration in the choice of a suitable model. 
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2. Extent of the Near-Fleld Region 


In determining the extent of the near-*fleld region, two aspects of the 
close-approach problem are of Interest: I) How close can a discrete vortex 

approach a discretized vortex sheet before its lateral position relative to the 
lattice control points distorts the vorticlty solution significantly? 2) How 
close can a discretized vortex sheet be approached before the calculated 
velocity distribution shows significant errors? The two aspects are clearly 
related and would be expected to yield similar near-field criteria, but it is 
instructive to consider them separat^^y, bearing In mind that both situations 
can occur In multiple-component calculations with force-free wakes. 

For this investigation errors greater than 0.5Z will be regarded as 
significant, l.e., this figure will define the edge of the near-fleld region 
here. It is well below current experimental accuracy for velocity s^asure- 
ments. 

2.1 Vorticlty Solution 

To examine aspect (1) of the close-approach problem, consider an infinite 
discrete vortex situated above and parallel to an infinite plane surface. The 
flow field is two-dimensional, and the vorticlty distribution induced on the 
surface is: 

Y(y) - -fZ/x[(y - Y)2 + Z^3 (1) 

where y is measured along the surface (Fig. 1), (Y, Z) is the vortex location, 
and r is the vortex strength. 

The region -1.0 <y 1.0 is now discretized using "N" equally spaced 
vortices. Control points are positioned midway between the discrete vortices 
(Fig. 1) following standard vortex- lattice procedure. The boundary condition 
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specified at the control points is Chet the vortices in the surfece induce 
the sene nomel velocity as that induced by the replaced part of the vortex 
sheet: 


M 


— 

2ir / ^ 


i-1 


ycj - yvi 




1 - 1» 2,.., N-1 


( 2 ) 


where y^ ■ jd - 1,0 are the control points; y^^ • (i - 0.5)d - 1.0 are tlw 
vortex points; A ■ 2.0/N is the vortex spacing; aad Ti are the (unknown) 


vortex strengths; 
and 

W(y) - - 


fz (i ( Fd .. - 

2»^t(y - If)* + 2^] ([a + Y)* + z*J\‘ ) 


The Nth equation concerns the conservation of circulation, i.e., the strengths 
of the surface vortices must add up to the circulation of the replaced part 
of the vortex sheet: 


i«l ^ 

Using Eqs. (2) and (3) we can solve for the vortex strengths and 
obtain the "discretized" vorticlty solution: 


(3) 


Yi - Fi/d 


In order to keep the edges of the discretized region remote from the 
region of interest (in terms of vortex spacing), a large m^er of discrete 
vortices (N ■ 40) were used. Solutions were obtained for a range of free 
vortex locations (Y, Z); at various heights above the sheet the onset vortex 
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was travaraad laterally from a atatlon above a control point to a station 
above a vortex point. 

The calculated vorticity distributions for three heights, Z * d, 0.5A, 
and 0.2A, are shown in Fig. 2 (centered on the onset vortex lateral position) 
together with the exact solution from Eq. (1). Appreciable deviations in 
the solution occur only below Z ■ A. Above that height the vortex position 
relative to the lattice arrangement is not critical, e.g., at Z • 2A the max- 
imuQ error in vorticity in the peak region is of the order of 0.(K)1Z, while at 
Z ■ A it is only 0.37%. With the onset vortex at Z ■ 0.5A the error near 
the vorticity peak varies from -f-IOZ to -10% as the vortex moves fr<» a point 
above a control point to one above a vortex point. When the vortex is at 
either end of the traverse, the vorticity errors decrease with distance from 
the peak vorticity position and are synm^trical about that point (Fig. 2b). 

For intermediate locations of the onset vortex, however, the errors are 
appreciable over a wide region, and they tend to be antlsymmetrical about 
the peak vorticity position. For very small heights, where the solution is 
tending towards a concentrated vortex, even larger antisymmetric errors 
appear in the solution unless the vortex lies above a control point or above 
a vortex point (Fig. 2c) . 

An interesting observation is that the solution with the vortex above a 
control point has a positive error, while that for the vortex above a vortex 
point has a negative error of about the same magnitude. Combining the two 
solutions at the middle station gives a very small error. It might be 
possible to make use of this fact in situations where a single vortex is very 
close to a discretized vortex sheet. 
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2.2 Velocity Dietrlbutlona 

To evaluate the near-fleld diatance for the velocity calculation aeprot* 
coneider a flat, two-dlaenalonal vortex aheet %rlth a parabolic vortlclty 
dletrlbutlon: 

Y<y) ■ y(l “ y) ; O < y < l 
The exact Induced velocity conqionenta are: 

V(y, z) - j[y(y - l) - Z^] T(y, z) + 2z(0.5 - y)L(y, z) + 

W(y, z) - |2a(0.5 - y) T(y, z) - [y(y - 1) - z^]L(y, z) + y - 0.5j^2ir 

where 

2 2 

LCy, z) • In (y » y^O, y^llfz-0 

and 

T(y, z) “ tan“^ tan”^ ~ — ) 

and V la the velocity component parallel to tlw vortex aheet; W la the normal 
component, l.e.. In the direction of z. 

Now dlacretize the vortex aheet ualng ”N" equally apaced vortlcea at 
locatlona: 

y^^ - (1 - 0.5)A ; 1-1,2 M 

where L - 1.0/N la the vortex apaclng. 

Integrating the vortlclty over each vortex region yielda the vortex 
atrengtha : 

ri - yvid.O - yvi)A - (aV12); 1-1, 2,...,N 

The velocity componenta Induced by the dlacretizatlon are: 

A h* 

N 

Wd(7. ^ 2 

1-1 
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Using forty disctsts iwrtlcss, velocity distributions were cslculstsd 
over s region betiMsn aldpoints near the quarter position on the sagaant 
(Fig. 3a). Error contours are shown in Fig. 3b. The discretisation gives 
negligible errors for both com>on«its of velocity in the region beyond Id froa 
the sheet. In effect the *'iwles" in the representation are not sensed tmtil 
we enter the Id region. 

Inside the Id region the errors Increase rapidly except along special 
lines of approach to the sheet. For the nomal velocity coaponent, the zero- 
error lines follow approxlaately the nomals to the sheet at the points aidway 
between the vortices and also at the vortices. (Deviations froa the normal 
lines occur because of the gradient in vorticity across the region.) Both 
sets of positions on the surface are used in the extended vortex- lattice 
method.^ l.e., the midpoints ar( used as control points, as in the standard 
vortex-lattice theory, and the vortex points ere used when applying the 
Kutta-Joukowski law for local forces and also when performing the tralling- 
vortex roll-up calculations. The zero-error lines for the tangential velocity 
component are less well known; these lines enter the near-fleld region above 
the quarter and three-quarter positions between the vortices and approach the 
vortex locations along approximately elliptical paths. All the zero-error 
paths are situated on extreme "precipices'* in the error contour map; small 
deviations from the paths result in large errors and lead to the near-fleld 
problems . 


3. Near-Field Models 

The previous section Indicates that errors arising from the discretiza- 
tion of a vortex sheet become appreciable only within the Id region. Clearly, 
if we wished to calculate velocities very close to the discretized vortex 
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sheet, ve could siaqply decreeee the Mice of A by increeslns the iiiad>er of 
vortices; however, for th*!:ee-diaenslofiel piobleBS the cowputing tiae could 
then becoM prohibitive. An eltemetlve solution is to epply e neer>field 
treetasnt to the wrtices. This trestswnt irauld be epplled only to those 
vortices thet ere within e specified "neer-field redlus” (e.g.. It) ir<m the 
point fdiere the velocity is being cslculeted. A ntn^er of neer-fleld aodels 
were considered in this work in two groups, i.e., "core-«odels'' end "spreed- 
■odels." ^e distinguish tog feeture betwe«i the groups is thet wherees core 
nodels applied without reference to neighboring vortex positions, spread 
fl»dels rec^ulre that Inforaation to generate connecting "vortex sheets." Cne 
nodel fr<» each group is described in this section. These lead to the sub- 
vortex tecltnlque, which is described separately in Sec. 4. 

3.1 Core Models 

A core model offers the simplest near-field treatment which remanws the 
singular behavior of the velocity field. In such a model the velocity induced 
by the vortex is factored locally so as to remain bounded at the vortex center. 
The Rankine Vortex and Laiid>'s viscous vortex are well known exa^les, but 
there are other possible forms. Core models have been used to the past to 
smooth the motions of vortices used to two-dimensional roll-up calculations 
(e.g.. Refs. 6 and 7). 

Several core models were tested using the discretized parabolic vortlclty 
sheet from Sec. 2.2,- but none were found sat le factory for both components of 
velocity. For example, they fall to restore the tengential com>onent of 
velocity near the vortex sheet. This can be seen to Fig. 3c (1) which shows 
the error contours for a Rankine vortex nodel with a core diameter of A. 
Although the tangential velocity errors appear slightly worse than for the 
utsMdifled vortex [compare Figs. 3c (1) and 3b], the normal component errors 
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At€ ^n>v«d on th« within th« corn. But thn nrror Itvnln are ■till 

■IffilflCMiti Md th« Baro*«rror lines no longer approach the vortex points. 

3.2 Spread Models 

In spread aodels the vortex Itself is aodlfiedi mad its strength is 
effective!/ distributed along a line representing the local position of the 
vortex sheet. Conpared with core aodels, spread nodela would be sore difficult 
to apply to the three-dinensional case and would be expected to consuoe nore 
cooputing time. 

There are a nuaber of ways of spreading the strengths of the wrtices; 
one based on a linear distribution of vortlcity proved very effectliw. In 
this aodel the >wrticity associated with a vortex is distributed in a tri- 
angular fashion on the two straight segn^nts joining the vortex to its 
innediate neighbors. When several neighboring vortices are treated in this 
way, the local effect is that of a piecewise, linear distribution of 
vortlcity— a model already known to work well in two-dimensional aerofoil 
theory (e.g. , Ref. 5). 

When evaluating this model on the discretized parabolic vortlcity distri- 
bution from Sec. 2.2, it was found necessary to extend the near- field radius 
of the vortices from Id to 4d to achieve the required accuracy. (The larger 
radius is required only for the lateral distance along the sheet; it is not 
required for the normal distance from the sheet, which is still Id.) The 
reason for the extended radius will be discussed in Sec. 4.2. 

With the larger near-field radius, the errors for this problem are less 
than 0.5Z everywhere in the flow field, including the surface of the vortex 
sheet. However, the three-dimensional form of this model can be cuiri»ersome 
to apply in a general method. A discretized form of this model was therefore 
considered, and this formed the basis of the subvortex technique. 
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4. Sttbvortcx TtchniqiM 


In thft fubvortex t«ctmiqu«, th« sprcnding of • noir>fi«ld vortox Is 
schisvsd by splitting it into s nus^r of sasll vorticss which sr« plsc^ 
bstwssn ths basic vortsx positions. Ths tsclmiqus offers s ntabsr of sdv«i- 
tsgss over ths linssr vorticity aodel uid has savaral novel faaturas* 

4.1 Features of the Technique 

(a) Subvortex positions 

The subvortices are distributed evenly along the vortex sheet joining the 
vortex to its two innediate neighbors. Unlike the linear vorticity aodel 
considered in Sec. 3.2, the joining sheet is not necessarily a straight line; 
the subvortices can be placed on an interpolated curve passing throuj^ the 
basic vortices, and so a better representation of curved vortex sheets is 
possible. Half intervals separate the basic vortices froa the nearest sub- 
vortices (Fig. 4a), and so the basic vortex positions becosM aidpolnts in the 
subvortex systea. This feature i^>roves the accuracy of the calculated 
velocity at the basic vortices (see Sec. 4.3). 

(b) Mtari>er of subvortices 

The number of suhvortlces used is such that the point iHiere the velocity 
is being calculated caimot "see the holes" in the discretised vortex slMet, 
l.e., the point is kept just outside the new local Id region of the aubvortex 
system. Figure 4b shows how this works using ths following expression for the 
nuaber of eubvortices on one side of the basic vortex: 

NSV - integer-part-of (1 + d/H) (4) 

where H is the noraal dlstuice of the poinc froa the aegaant. Use of this 
expression keeps the nuaber of subvortices to s sdnisaai and helps to keep 
cosq>utlng costs down. Uhen applied to the vortex-lattice methods, the add- 
points between the vortices (l.e., the cmtrol points) should rMuiin 
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aldpointf in th« tubvortnx syntmii; N8V wu»t th«i b« tvn . «• shown 

dotted in Fig. 4b* 

A mutimm llait is plscsd on ths msi^r of subwrticss to svoid s 

runswsy condition when ths height H spproschss ssro. This lialt controls ths 
closest spprosch thst csn be asds to the vortex sheet before the new local 
Id region of the subvortices is entered. It csn therefore be used to control 
the '*eccurscy” of the cslculstion in s trsde*off with coeluting tisM, i.e., 
by incrsssittg the liait the error regira lenild decrease in sice, but the 
coBg>uting tiM would increase, and vice-i^rss. 

(c) Subvortex strengttw 

The subvortices aust have s coabined strength equal to thst of the 
sssociste^. basic vortex. In the technique as used here their strengths 
vary linearly with distance froa the basic vortex position. With equal 
spacing, therefore, the 1th subvortex has strength 

r. - T tiLo T - 0.5) ; for i • 1, 2,..., I«V on each side 

1 Qf basic vortex. Fig. 4a. 

where f is the strength of the basic vortex. When several neighboring basic 

wrtices are treated in this way, the local effect approaches that of a 

piecewise linear vorticlty distribution. Clearly, higher order distributiwis 

could be used but would involve wore than one basic vortex Interval on each 

side. 

(d) Subvortex cores 

Because of the practical lialt placed on the nua^r of subvortices, 

it is possible for a point irtiere the velocity is calculated to fall within 
the (reduced) Id region of the subvortex syst«s. Because of this each 
subvortex has been aodifled with a Rankine-Vortex core (Sec. 3.1) of diaawter 
approxiaately equal to the distMce between the subvortices. Tills saears out 
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the dl8c<mtiiM'lty associated with the sheet but mily ovbt the new, dlBinished, 
Mar-field region. 

An alternative, but wore ccwplicated possibility a subvortex near- 
field region is entered is to apply a special “sheared core." For this the 
velocity point is physically moved the snail distance to the nearest midpoint 
between subvortices when evaluating the local contribution. The usual velo- 
city calculation then provides the local mean velocity, to which is added the 
local vortlcity contribution. The latter is half the local vorticity value 
(evaluated by dividing the nearest subvortex strength by the subvortex spacing) 
and is directed along the local tangent on one side of the vortex sl^t and 
in the opposite direction on the other side. A similar sheared core was one 
of the models tried earlier (Sec. 3. 1) but was not satisfactory for the broad 
intervals of the basic vortices; hoMver, it should lq>rove the subvortex 
technique. 

(e) Extension to three dlieensions 

The subvortex tectmlque can be adapted to the segmented vortices used in 
vortex-lattice-based methods. The location of the ends of each subvortex 
segment is a sitq>le extension of the two-dimensionax method. 

4.2 Near-Field Radius 

As in the case of the linear vortlcity model (Sec. 3.2), when testing the 
subvortex technique, the near- field radius had to be Increased beyond lA to 
obtain the required accuracy. The reason for this extension is that the 
Induced velocity from the "distributed" model does not match that from the 
basic vortex until some distance away. For the subvortex model this result is 
fairly insensitive to the nua^er of subvortices (see Fig. 5). Thus, when one 
vortex has the near-field modification applied, several neighboring vortices 
must also be modified, otherwise small "jumps" occur in the calculated velocity 
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distribution as we pass over the near-field boundary of each vortex. These 
jiaaps can be aade as sautll as we please by increasing the near-field radius. 
For the subvortex technique the jiuq>s are effectively elisinated whm using a 
near- field radius of 5A. 

Although the near- field radius of 5d is strictly required only in the 
direction parallel to the vortex sheet (the nomal distance is still lA), a 
true radius of SA is now used. This can result in superfluous calculations if 
the point actually lies beyond lA from the sheet, but it does provide more 
reliable conditions for highly cun^d vortex sheets. The possible penalty in 
coflq>uting time is alleviated by the fact that the tangential distance between 
the point and the segment is no longer required and that only a small ntu^r 
of subvortices are used for the outer region, l.e., the region lA to 5A 
(Fig. 4b). 

4.3 Error Contours 

The technique was tested on the discretized parabolic vortlcity distri- 
bution considered In Sec. 2.2. The error contours (Fig. .3c(ll)) appear not 
quite as good as those for the linear vortlcity model (Sec. 3.2); a small 
region with significant error is indicated very close to the vortex sheet 
where the approach is closer than the subvortex spacing and is associated 
with the maximum limit placed on the number of subvortices (Sec. 4.1b). In 
these calculations NSV^^^ was 10. 

The normal component of velocity calculated at the vortex locations has 
always been slightly less accurate than that calculated at points midway 
between the vorUices. (The \»ortex points are effectively midpoints in a 
coarser discretization.) For the present discretized parabolic vortlcity 
distribution, the error at the vortices in the region considered (sec Fig. 3a) 
is 2.8% compared with 0.03% at the midpoints. With the subvortex technique 
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applied, the error at the vorticea decreaaea to 0.2%; this reduction la helped 
by the fact that the basic vortex locations bectxRe midpoints In the subvortex 
system (Sec. 4.1a). 


5. Two-Dimensional Application 

For a more general situation the subvortex technique vaa applied to the 
two-dimensional, time-dependent, roll-up calculation of the vortex sheet shed 
from a wing with elliptical load distribution: 

Y(y) • 2y/ yi - y^ ; -1 ^ y < 1 

This problea has been considered by a number of authors, e.g. , Refs. 6-8 
(most recently). The vortex sheet Is first discretized then the motion of the 
discrete vortices under their mutual Interactions is followed over a number 
of small time steps fit. After the 1th time step a point on the trajectory of 
a vortex Is calculated here using the expression 

li + 1 ■ 

where Is a mean velocity calculated at a midpoint in the time interval. 

This point is obtained by extrapolating the conditions at the two previous time 
steps. Before calculating all the vortices are moved to their midpoint 
positions. Unlike most previous roll-up calculations, the Influence of the 
vortex on Itself is not Ignored here; because the vortex point Is a midpoint 
In the subvortex system (Sec. 4. la) , the vortex can now Influence Itself if 
the vortex sheet is curved or stretched locally. 

In the present calculation ten equally spaced vortices are used on one 
half of the sheet, and the time step fit is 0.075. (This is smaller than the 
shortest orbital period of the vortex pairs In the system (see Ref. 8), that 
value being 0.2 here.) The vortex positions are shown in Fig. 6a after 
twenty time steps, l.e., after a tliim of 1.5. The curve drawn through the 


- 14 - 



vortex locations Is the interpolated curve on which the subvert icea are 
placed. With so few vortices in the roll-up region. Interpolation there is 
difficult; the scheme used here is based on a lliwar combination of two 
circular arcs over each Interval, but even this breaks down (i.e., crosses 
Itself) after a few more time steps beyond that shown. (When developing the 
three-dimensional version of the model, such problems will be alleviated by the 
use of a vortex^amalgamation technique, such as that used by Moore. 

Hsvlng obtained a representative distribution of discrete vortices, a 
velocity scan was taken through the roll-up region at y 0.8. The •'wo 
velocity components are shown in Fig. 6b, c with and without the subvortex 
technique being used. The basic discretization (with small Ranklne-vortex 
cores applied) shows the characteristic velocity deviations near the discrete 
vortices and Ignores the presence of connecting vortex sheets. Qualitatively, 
the subvortex technique shows a plausible modeling of the discontinuities 
associated with the (Inviscld) vortex sheet. By Increasing the limit on the 
number of subvortices (from the ten used here), the representation of 

the discontinuities would be sharpened. Quantitative comparisons are still 
required for a general case, but these will be pursued after the technique 
has been Incorporated In a three-dimensional method. 

The end vortex near the center of the roll-up region requires further 
treatment, possibly In the form of a viscous core; such a treatment coupled 
with the analgamatlon technique mentioned earlier should Improve the cal- 
culations near a roll-up region. 


6. Conclusions 

Discretization of a vortex sheet Introduces significant velocity errors 
only within a distance from the sheet equal to the vortex spacing In the 
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lattice. Ck>n nodele applied to the vorticea help to limit the aiee of errore 
but do not reduce them to a aatisfactory level when the field of interest 
approaches close to the vortex sheet. The region where significant errors 
occur c«n be reduced to a small region of controllable width close to the 
vortex sheet by the use of the near^field model in which a discrete vortex 
splits into an increasing nusd>er of subvortices as it is approached. The two- 
dimensional probleiM considered here dononstrate that the subvortex technlqtw 
can achieve the objective of this investigation, viz., to provide accurate 
velocities anywhere in the flow>fleld irrespective of the proximity of dis- 
cretized vortex sheets. (Clearly, quantitative ctmiparlsons with eiq>erim«ital 
results are still required, but these will be pursued after the model has 
been Incorporated in a three-dimensional method.) 

In sumaary, the subvortex technique promises to enhance the versatility 
of vortex-lattice-based methods by providing the effect of a much finer 
discretization. The technique is efficient because the ntmd>er of vortices in 
the lattice is effectively Increased but only where and when needed and by an 
amount just sufficient to prevent the ’’holes" In the lattice from being "seen." 
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(a) Vortex height, Z • L, 

Figure 2.- Vortlcity Induced on a plane surface by a vortex — effect of vortex 

location on dlacretlzed solution. 



(b) Vortex height, Z • 0.5A 


VORTEX LOCATION 

Z/A=0.2 



(c) Vortex height, Z > 0.2A 





Figuni 3.- Velocity error celculetlons 




(11) Subvortcx t«chnlquc; ncxr-f d radius * 5d» WSV^^ ■ 10. 


(e) Error contours for tvo nssr-finld ■odsls. 
Flfuro 3.» Concluded. 
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(b) Niisd)er of subvortices as a function of height above the segment. 


Figure 



5.- Subvortex technique. Difference between velocities from subvortices 

and from the basic vortex. 



(a) Calculated vortex positions after 20 time steps (t - 1.5). 

Figure 6.- Single application of the subvortex technique to the vortex sheet 

roll-up calculation In two dimensions. 
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(b) Calculated upwaeh distrlbutlcm 


SIDEWASH VELOCITY, 



(c) Calculated sidewaah distribution. 
Figure 6.- Concluded. 


